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SUMMARY 


Modern  numerical  methods  for  the  optimization  of  large  dis¬ 
cretized  systems  are  now  we i L  developed  and  highly  efficient  in 
the  case  of  thin  vj ailed  elastic  structures  mcaeled  by  finite 
elements.  However,  this  is  not  yet  true  for  structures  whose 
components  are  subject  simultaneously  to  bending  and  extension 
loads.  In  this  paper,  the  idea  of  Generalized  Optimality  Cri¬ 
terion  (GOC),  set  forth  in  previous  final  scientific  reports 
for  bar,  membrane  and  pure  bending  elements,  is  extended  to  deal 
with  general  beam  and  flat  shell  elements.  The  modifications 
brought  to  the  GOC  result  in  explicit  approximations  for  the 
behavior  constraints  that  are  still  correct  up  to  the  first  order, 
but  that  exhibit  a  more  comp,  lex  algebraic  form.  Indeed  these 
explicit  expressions  are  no  longer  merely  linear  in  the  recipro  .il 
design  variables.  However  they  continue  to  be  addi lively  sepa¬ 
rable  and  therefore,  dual  methods  remain  fully  applicable,  just 
as  in  the  original  statement  of  the  GOC  approach.  Numerical 
examples  will  be  offered  to  demonstrate  the  efficiency  of  the 
method  presented. 
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INTRODUCTION 


The  optimum  design  of  any  significant  structure  is  the  result 
of  a  delicate  compromise  between  many  complex  factors.  Some  art- 
rational  and  can  be  quantified,  such  as  the  strength  of  t fit- 
structure,  its  natural  frequencies,  its  weight,  ...  Some  are 
just  as  rational  but  are  difficult  to  quantify,  such  as  the  ex¬ 
perience  of  the  designer  in  a  given  technology.  Some  others  are 
much  less  rational,  like  styling,  but  are  just  as  important  for 
the  final  goal  of  the  process,  which  is  the  marketing.  Naturally 
a  good  designer  considers  structural  optimization  as  a  technique 
that  should  take  into  account  all  possible  aspects  of  the  design. 
Consequently  the  designers  are  often  reluctant  to  the  concepts 
of  structural  optimization  developed  in  connection  with  finite 
element  programs. 

However  a  more  detailed  examination  of  the  design  process 
leads  to  isolate  a  phase  that  appears  frequently,  during  which 
the  shape  of  the  structure  is  more  or  less  frozen  and  the  problem 
is  limited  to  giving  adequate  dimensions  to  the  various  members. 
Such  a  situation  is  often  encountered  in  the  aerospace,  naval 
or  automobile  industries,  where  the  external  shapes  are,  to  a 
large  extent,  dictated  by  aero-  or  hydrodynamic  considerations, 
or  by  styling,  while  internal  forms  are  often  determinated  by 
various  other  non  structural  requirements.  If  the  ultimate  goal 
of  the  designer  can  be  identified  as  corresponding  to  the  minimi¬ 
zation  of  an  explicit  function  of  the  member  sizes,  and  if  the 
limitations  on  the  design  can  be  defined  as,  eventually  implicit, 
functions  of  the  member  sizes  too,  such  as  displacements,  stresses, 
eigenfrequencies,  etc...,  then  the  problem  is  tractable  by  auto¬ 
matic  algorithms.  They  allow  the  designer  to  speed  up  signifi¬ 
cantly  this  part  of  the  design  process  and  to  explore  more  sys¬ 
tematically  the  various  feasible  designs. 

The  optimization  problems,  which  in  fact  should  be  called 
automatic  sizing  problems,  are  especially  crucial  when  complex 
structural  forms  are  involved,  when  flexural  forces  cannot  be 


2  . 

neglected,  and  when  composite  materials  such  as  reinforced 
resins  are  employed.  In  these  cases  it  becomes  difficult,  if 
not  impossible,  for  the  designer  to  have  an  intuitive  under¬ 
standing  of  the  structural  mechanics  that  is  sufficient  to  leu! 
to  op t i mal  s i z i ng  of  the  various  members.  Furthermore,  the  de¬ 
signer  is  most  of  the  time  unable  to  take  into  account  global 
constraints  in  the  structure,  like  global  flexibility,  restric¬ 
tion  on  displacements,  frequencies  of  vibration,  global  buckling 
modes,  etc...  It  is  only  possible  to  verify  a  posteriori  that 
such  constraints  are  satisfied.  Again  these  global  constraints 
become  more  important  in  the  context  of  highly,  indeterminate 
structures.  In  the  aerospace  industry,  the  necessity  of  desi¬ 
gning  high  performance  structures  has  motivated  significant  re¬ 
search  efforts  to  derive  algorithms  permitting  a  rapid  and  sys¬ 
tematic  exploration  of  the  design  space  to  determine  the  optimum 
material  utilization. 

It  is  worth  pointing  out  that  optimization  methods  should 
be  considered  as  especially  useful  in  the  preliminary  design 
phase.  Using  them  when  the  design  is  practically  frozen,  with 
the  hope  of  an  ultimate  improvement,  is  often  disappointing. 

This  i  due  to  the  fact  that  the  optimization  of  a  detailed  design 
implies  the  formulation  of  a  large  number  of  constraints,  some 
of  which  are  not  easily  quantified.  At  the  preliminary  design 
stage,  however,  the  constraints  are  usually  more  global  and 
therefore  more  easily  handled  by  the  available  formulations. 

The  structural  optimization  problem  considered  in  this  report 
consists  of  the  weight  minimization  of  a  finite  element  model 
with  fixed  geometry  and  material  properties.  The  design  variables 
are  taken  as  the  transverse  sizes  of  the  structural  members, 
namely,  the  cross-sectional  areas  of  bar  and  beam  elements  and 
the  thicknesses  of  membrane,  plate  and  flat  shell  elements. 

The  mathematical  programming  problem  to  be  solved  has  the  fol¬ 
lowing  form  : 
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The  aVs  denote  the  n  design  variables.  They  correspond  to 
member  sizes  of  either  individual  finite  elements,  or,  if 
design  variable  linking  is  used,  of  groups  of  finite  elements. 
The  structural  weight  W  is  a  linear  objective  function,  because 
the  £  ^  '  s  are  constant  coefficients  representing  the  specific 
weight  times  either  the  element  length  (bars  and  beams)  or  the 
element  area  (membranes,  shearpanels,  plates  and  flat  shells). 
The  inequalities  (2)  are  the  behavior  constraints,  which  impose 
limitations  on  quantities  describing  the  structural  response, 
for  example,  the  stresses  and  displacements  under  multiple 
loading  cases,  the  natural  frequencies,  the  buckling  loads, 
etc...  The  design  variables  are  also  subjected  to  the  side 
constraints  (3),  where  a^  and  a^  are  lower  and  upper  limits  that 
reflect  fabrication  and  analysis  validity  considerations. 

Standard  minimization  techniques  can  be  applied  to  the  non¬ 
linear  programming  problem  (1-3).  However  this  problem  exhibits 
some  characteristics  that  make  it  complicated  when  practical 
structural  design  applications  are  considered.  The  essential 
difficulty  arises  from  the  implicit  nature  of  the  behavior 
constraints  (2),  in  that  their  precise  numerical  evaluation  for 
each  particular  design  requires  a  complete  finite  element  ana¬ 
lysis.  Since  the  solution  scheme  is  iterative,  it  involves  a 
large  number  of  structural  reanalyses.  Therefore  the  computa¬ 
tional  cost  often  becomes  prohibitive  when  large  structural 
systems  are  dealt  with.  However  a  powerful  design  procedure 
has  now  emerged,  which  consists  in  replacing  the  initial  problem 
with  a  sequence  of  simple  explicit  problems.  In  the  next  sec¬ 
tion  this  approach  will  be  briefly  reviewed  by  restricting  the 
formulation  to  thin  walled  structures  idealized  by  bar/membrane 
elements.  A  mucli  more  detailed  presentation  can  bo  found  in  a 
previous  report  II].  It  will  be  shown  that  the  behavior  cons¬ 
traints  can  be  approximated  either  by  using  virtual  load  consi¬ 
derations  (optimality  criteria  approach)  or  by  using  first  order 
Taylor  series  expansion  with  respect  to  the  reciprocal  design 
variables  (mathematical  programming  approach).  Applying  a  dual 
solution  scheme  to  each  explicit  problem  generated  in  sequence 
naturally  introduces  the  concept  of  a  generalized  optimality 


criterion. 


S ub s eq ue n t 1 y s e c t i on s  3  and  4  wilL  be  concerned  with  struc¬ 
tural  models  that  are  capable  of  carrying  flexural  forces.  For 
beams  and  plates  in  pure  bending,  adequate  intermediate  varia¬ 
bles  can  be  selected,  in  terms  of  which  high  quality  explicit 
approximations  for  the  behavior  constraints  can  still  be  gene¬ 
rated  by  linearization.  The  idea  of  generalized  optimality 
criterion  remains  fully  valid  and  it  keeps  its  interpretation 
in  terms  of  energy  densities  in  the  structural  members.  Section 
3  is  a  summary  of  results  presented  in  a  previous  report  (ref. 

[2  j,  section  6). 

The  essential  problem  is  addressed  in  section  4.  It  consists 
of  the  establishment  of  the  generalized  optimality  criterion 
approach  in  the  general  case  where  the  structural  members  work 
both  in  extension  and  flexion  (beam  and  flat  shell  elements). 

For  displacement  constraints,  it  is  no  longer  possible  to  select 
a  suitable  intermediate  variable  for  the  linearization  process, 
however,  the  virtual  load  procedure  permits  obtention  of  high 
quality,  first  order  explicit  approximations  of  the  behavior 
constraints.  It  will  also  be  shown  how  to  proceed  for  stress, 
frequency  and  buckling  constraints. 

The  explicit  approximations  still  exhibit  a  separable  alge¬ 
braic  form  and  therefore  dual  methods  remain  applicable.  However 
some  difficulties  might  happen  due  to  the  lack  of  convexity  in 
the  explicit  subproblera.  Section  3  will  describe  how  to  solve 
the  explicit  subproblems  by  using  dual  methods.  Finally  some 
applications  will  be  offered  in  section  6  to  illustrate  the  power 
and  generality  of  the  approach  presented. 


2  .  GENERALIZED  OPTIMALITY  CRITERION  FOR  THIN  WALLED  STRUCTURES 

This  section  summarizes  some  results  obtained  in  a  previous 
work  (  1  J  for  structural  models  made  up  of  bar,  membrane  and 
shear  panel  elements,  which  are  quite  adequate  for  idealizing, 
thin  walled  structures  subjected  mainly  to  extension  loading. 
For  this  class  of  finite  element  models,  the  structural  stif¬ 
fness  matrix  exhibits  a  linear  form  in  the  design  variables  : 


K  = 


K  . 
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i  =  I 


a  . 
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( '•  ) 


where  K^,  a  matrix  of  constants,  represents  the  stiffness  matri 


of  the  ith  element  when  a  ^  =  i 


For  simplicity,  the  following, 
discussion  is  restricted  to  problems  involving  constraints  on 
static  stresses  and  displacements,  in  which  case  the  behavior 
constraints  (2)  can  be  written 
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u  . 
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u . (a)  >  0 
J 


(  ')  ) 


where  u.  denotes  an  upper  bound  to  a  response  quantity  Ujfa) 
(stress,  nodal  displacement,  relative  displacement). 


Most  of  the  optimality  criteria  approaches  (e.g.  |  3  1 ),  as 
well  as  the  generalized  optimality  criterion  (GOC)  set  forth  in 
Ret.  [11,  use  the  virtual  load  technique  to  generate  explicit 
approximations  of  the  stress  and  displacement  constraints.  In¬ 
troducing  a  virtual  load  vector  conjugated  to  the  response  quan¬ 
tity  Uj  (unit  load  for  a  nodal  displacement),  it  follows  that 
u j  can  be  expressed  as  the  sum  of  the  contributions  of  each  fi¬ 
nite  element  : 


T  n  Ci  j 

uj  =  q  K  q.  =  .£,  - i 
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with 
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Ln  these  expressions  q  and  are  respectively  the  real  and 

virtual  displacement  vectors  and  K  .  is  the  element  s  t  i  I'  fine,;; 

matrix  appear!  hr  in  (4).  It  can  be  seen  from  (7)  that  the  •  i>e!- 

ficients  e .  .  are  related  to  the  virtual  strain  one  rev  densities 
i  J 

in  the  structural  members.  The  c..'s  are  constant  coefficients 

i  J 

in  the  case  of  a  statically  determinate  structure,  set  that  (<>) 
represents  then  the  exact  explicit  form  of  the  response  quan¬ 
tity  Uj.  In  the  case  of  a  statically  indeterminate  structure, 
the  c^j's  depend  implicitly  on  the  design  VouiabJcJ,  because 
structural  redundancy  produces  redistribution  of  the  internal 
forces  when  the  member  sizes  are  modified.  Therefore  the  lol- 
1  owing  explicit  constraints 


'v  -  n  i  i 

h  .  ( a )  u  .  -  .  ,  - — J-  >  0  18) 

J  J  i=J  a. 

constitute  in  general  approximate  forms  of  the  original  cons¬ 
traints  (5).  As  shown  in  Ref.  !  1  ] ,  the  basic  idea  in  the  op¬ 
timality  criteria  approach  can  be  viewed  as  transforming  the 
initial  implicit  problem  into  a  sequence  of  explicit  subproblem 
Each  explicit  problem  results  from  replacing  the  behavior  cons¬ 
traints  (2)  by  their  approximate  forms  (8). 


On  the  other  hand,  the  mathematical  programming  approach  to 
structural  optimization,  after  a  period  of  unefficiency,  has 
finally  evolved  into  a  powerful  and  now  well  established  design 
procedure  which  is  also  based  upon  explicit  approximations  o I 
the  behavior  constraints  [4,  3,  6  ].  The  key  idea  is  to  linea¬ 
rize  the  behavior  constraints  with  respect  to  the  reciprocal 
design  variables 


Justification  for  this  change  of  variables  lies  in  the  fact 
that  the  constraint  surfaces  can  be  shown  to  be  very  shallow 
and  close  to  planes  in  the  reciprocal  design  variable  space. 
Therefore  the  linearized  forms  of  the  constraints  are  usually 
high  quality  approximations.  They  are  obtained  by  using  first 
order  Taylor  series  expansion  in  terms  of  the  reciprocal  va¬ 
riables  x.  : 
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n  ‘ u  •  ° 
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where  the  up  pcrsrript  denotes  quantities  evaluated  it  the 
actual  design  point  x°,  wliere  the  structural  analysis  is  per¬ 
formed.  Note  that  tli  e  finite  element  analysis  capability  must 
include  auxiliary  sensitivity  analyses  for  evaluating  the  lirsi 
partial  derivatives  of  the  response  quantities.  Most  often  the 
well  known  pseudo-loads  technique  is  employed  |  7  ] . 

It  has  been  shown  in  a  previous  report  |  I  |  that  the  explicit 
approximations  of  the  behavior  constraints  used  in  both  the 
optimality  criteria  and  mathematical  programming  approaches  ( p.  q  b 
8  and  10,  respectively)  are  identical.  Indeed  the  virtual  si  rain 
energy  densities  employed  in  the  optimality  criteria  appro¬ 

aches  are  nothing  else  that  the  gradients  of  the  response  quan¬ 
tities  with  respect  to  the  reciprocal  variables  : 
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Furthermore  the  definition  of  the  c..'s  following  from  the  vir- 

i  J 

t u a  1  load  technique  (see  E  q  .  7)  clearly  indicates  that 
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Therefore  (10)  can  bo  rewritten 
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which  is  equivalent  to  (8)  when  recast  in  terms  of  the  direct 
variables  a^.  It  is  thus  apparent  that  a  unified  approach  to 
structural  weight  minimization  of  finite  element  systems  has 
emerged,  which  consists  in  replacing  the  initial  problem  (1-3) 
with  a  sequence  of  explicit  approximate  -  or  linearized  -  pro¬ 
blems  of  the  following  form  : 
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In  order  to  compute  the  lagrangian  multipliers  satisfying  (20- 
21),  an  interesting  approach  is  to  resort  to  dual  methods,  which 
leads  to  maximize  the  lagrangian  function  considered  as  a  func¬ 
tion  of  the  lagrangian  multipliers  -  or  dual  variables  -  only 
(see  Refs  [I,  6,  8  ]  and  section  5).  Once  the  solution  of  the 

dual  problem  has  been  found,  the  corresponding  optimal  design 
variables  are  easily  computed  from  the  explicit  optimality  cri¬ 
teria  relations  (17-19).  Note  that  the  design  variables  can  be 
separated  into  a  group  of  active  (or  free)  variables  (see  Kq  .  17) 

and  a  group  of  passive  (or  fixed)  variables  (see  Eqs.  18  and  19). 
This  subdivision  into  active  and  passive  design  variable  groups 
is  classical  in  the  optimality  criteria  approaches  |  3,  9  ]  .  It 

corresponds  to  the  fact  that  the  dual  space  -  i.e.  the  space  of 
the  lagrangian  multipliers  r^  -  is  partitioned  into  several  re¬ 
gions  separated  by  planes  across  which  the  second  derivatives 
of  the  dual  function  are  discontinuous  [6,  8  ] . 


The  whole  process  of  combining  the  linearization  of  the 
behavior  constraints  with  respect  to  the  reciprocal  design 
variables  and  a  dual  solution  scheme  can  be  viewed  as  a  gene¬ 
ralization  of  the  optimality  criteria  approaches.  It  is  impor¬ 
tant  to  mention  that  this  basic  approach  of  converting,  the 
initial  problem  into  a  sequence  of  explicit  s  ub  p  r  <>b  1  eins  is  now 
widely  recognized  [6,  10  ]  and  it  is  routinely  employed  for 

large  scale  industrial  applications  [  11  ]  . 


'3  .  P  U  R  E  B  E  N  D  [  NG  ELEMENTS 

In  this  section  attention  is  focused  on  discretized  models 
made  up  of  pure  beam  and  plate  elements  subjected  to  flexural 
loads  only  (for  more  details,  see  section  6  of  Ref.  |2  J). 

The  stiffness  matrix  of  such  a  bending  element  is  usually  not 
merely  proportional  to  its  cross-sectional  size  and  therefore 
the  optimization  strategy  reviewed  in  the  previous  section  must 
be  modified.  The  way  to  deal  with  a  beam  element  subjected  to 
uniaxial  bending  depends  upon  the  relationship  between  the  prin¬ 
cipal  moment  of  inertia  I  and  the  cross-sectional  area  a.  A 
wide  variety  of  situations  is  taken  into  consideration  by  adop¬ 
ting  the  following  relation  : 

I  =  c  aP  (22) 


where  c  is  a  constant  that  depends  only  on  the  shape  of  the 
beam  cross-section  and  p  is  a  positive  number. 


Most  of  the  time  p  is  taken  as  an  integer  number,  equal  to 
1,  2  or  3.  The  case  p  =  1  corresponds  to  thin  walled  beams, 
for  example,  sandwich  beams,  pipes  with  fixed  diameter  and  va¬ 
riable  thickness,  etc...  The  GOC  approach  of  section  2  remains 
then  fully  applicable,  since  the  stiffness  matrix  continues  to 
be  linear  in  the  design  variables.  The  case  p  =  2  is  that  of 
beams  with  uniformly  varying  cross-section.  The  shape  of  the 
cross-section  is  kept  unvariant  while  its  area  is  modified  during 
redesign  (dilatation  or  contraction).  Finally  the  case  p  =  3 
is  concerned  with  beams  having  full  cross-section  whose  height 
varies  while  other  sizes  are  fixed.  For  a  beam  subjected  to 
pure  bending,  the  flexural  rigidity  is  proportional  to  the  mo¬ 
ment  of  inertia  and  therefore,  in  a  finite  element  context,  the 
structural  stiffness  matrix  exhibits  the  following  explicit 
form  in  terms  of  the  cross-sectional  areas  : 
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where  each  matrix  K  ^  is  independent  of  the  design  variables 
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With  regard  to  plate  elements  subjected  to  pure  bending, 
two  cases  must  be  distinguished.  The  first  case  is  that  of 
sandwich  plates  with  constant  core  thickness.  The  sheet  thick¬ 
nesses  constitute  then  the  design  variables.  Consequent  ly  the 
stiffness  matrix  continues  to  depend  linearly  on  the  design 
variables  and  the  GOC  relations  (17-21)  remain  fully  applicable 
The  second  case  is  concerned  with  full  plates  with  variable 
thickness.  The  stiffness  is  then  proportional  to  the  cube  of 
the  thickness  and,  in  an  assembling  of  plate  elements,  the  re¬ 
lation  (22)  must  be  chosen  with  p  =  3. 

By  assuming  that  the  structural  discretization  is  made  up 
entirely  of  elements  of  the  same  type,  the  stiffness  matrix 
exhibits  the  form  (23),  where  p  takes  on  the  same  value  frn 
each  member.  In  these  circumstances,  the  GOC  can  be  derived 
just  as  in  the  case  of  thin  walled  structures,  by  adopting  a 
change  of  variables  tending  to  reduce  the  nonlinear  character 
of  the  constraints  : 


a . 
r 

The  next  step  is  to  linearize  the  constraints  with  respect  to 
the  new  variables  x^,  which  requires  gradient  evaluation  (see 
E  q  .  10).  Restricting  again  the  discussion  to  stress  and  dis¬ 
placement  constraints,  it  is  easily  shown  (see  Ref.  I  2  |)  that 
tfie  first  order  Taylor  series  expansion  (10)  reduces  to  the 
form  (13),  or,  when  written  in  terms  of  the  direct  variables  : 
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The  c^j  coefficients  can  be  interpreted  as  the  gradients  of  the 

response  quantities  with  respect  to  the  intermediate  variables 

x^  defined  in  (24),  but  they  can  also  be  related  to  the  virtual 

strain  energies  e. .  in  the  structural  members  : 
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where,  by  definition. 
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>  1  i  l 

In  this  connection  it  should  he  recognized  that  virtual  load 
considerations  could  directly  be  employed  to  derive  the  explicit 
approximations  (25),  instead  of  resorting  to  first  order  Taylor 
series  (see  next  section). 


The  optimality  criteria  equations  are  very  similar  to  (17- 
11),  provided  that  care  is  taken  of  the  exponent  p  appearing, 
in  (25).  The  basic  redesign  relations  for  the  active  design 
variables  must  read  as  follows  : 
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where  it  is  understood  that  the  dual  variables  r.  (i.e.  the 
lagrangian  multipliers)  must  satisfy  the  complementarity  con¬ 
ditions  (see  Eqs.  20  and  21).  A  physical  interpretation  of  the 
optimality  criterion  is  obtained  by  introducing  the  virtual 
strain  energy  densities  per  unit  weight  : 
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where  e..  is  defined  in  (27).  in  terms  of  the  <  .  .  ’  s  ,  the  o  p  t  i 
tj  i) 

mality  criterion  (28)  takes  the  "energetic"  form 
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In  the  special  case  where  only  one  displacement  constraint 
is  specified,  the  optimality  criterion  states  that  the  virtual 
strain  energy  density  must  be  the  same  in  each  element.  In 
this  simple  case,  it  is  possible  to  solve  analytically  the  ex¬ 
plicit  problem  and  to  derive  explicit  redesign  relations  in  terms 
of  known  quantities.  The  active  design  variables  can  be  shown 
to  be  given  by 
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while  the  remaining  passive  variables  are  fixed  to  an  upper  or 
a  lower  limit  (  uq  denotes  the  contribution  of  these  passive 
variables  to  the  displacement  constraint  u  <  u).  It  is  worth 
mentioning  that  (31)  is  well  suited  for  the  design  o  f  plates 
in  bending  with  a  single  displacement  constraint.  Since  then 
p  =  3  the  redesign  relation  (31)  involves  the  fourth  root  of  the 
coefficients  c  ^ ,  rather  than  the  third  root  as  employed  in 
Ref.  [12  1  on  an  intuitive  basis.  Note  also  that  by  taking 
p  =  1  in  (31),  conventional  redesign  relations  are  recovered, 

which  were  devised  for  trusses  [3  ]  ,  sandwich  beams  |  I  i  j  .etc... 

It  can  be  concluded  that  the  generalized  optimality  criteria 
approach  can  easily  be  extended  to  deal  with  pure  bending  ele¬ 
ments  by  defining  adequate  intermediate  variables.  High  quality 
explicit  approximations  of  the  behavior  constraints  can  s  t i  1  1  be 
generated  and  the  resulting  GOC  keeps  its  interpretation  in  terms 
of  energy  densities  ir  the  structural  members.  Seeking  the  de¬ 
sign  variables  that  satisfy  the  GOC  at  each  redesign  stage  can 
still  be  achieved  efficiently  by  resorting  to  dual  methods,  be¬ 
cause  the  explicit  approximate  problem  remains  separable  and 
strictly  convex  when  expressed  in  the  intermediate  design  vari¬ 
ables  x  ^  . 


4.  FLEX 1 ON-E  XTEN  S ION  ELEMENTS 


When  flexion  and  extension  loadings  act  simultaneously  witli 
comparable  intensity  at  the  element  level,  the  definition  (21) 
of  the  stiffness  matrix  can  no  longer  characterize  the  structu¬ 
ral  model  with  sufficient  accuracy.  To  help  fix  ideas,  consi¬ 
der  a  flat  shell  element  made  up  of  a  membrane  and  a  plate  sta¬ 
cked  together.  The  stiffness  matrix  of  such  a  Mat  shell  ele¬ 
ment  exhibits  the  form 


K.  . 
1 


♦  a3  K.(3) 
t  i 


(  32  ) 


where  and  K  f  3  ^  are  constant  matrices.  As  a  result,  in 

t  l 

the  GOC  approach,  if  the  constraints  are  linearized  with  respect 
to  the  reciprocal  design  variables  (9),  their  first  order  expli¬ 
cit  approximations,  given  by  expressions  similar  to  (8),  will 
be  of  high  quality  only  if  the  structural  members  behave  mainly 
in  extension.  On  the  other  hand,  if  the  bending  behavior  is 
dominant,  it  is  better  to  adopt  the  change  of  variables  (24), 
yielding  first  order  explicit  approximations  of  the  form  (2'i) 
(with  p  =  3).  As  a  matter  of  fact,  the  true  situation  is  usu¬ 
ally  a  combination  of  extension  and  bending.  In  a  practical 
structure,  some  members  work  mainly  in  extension,  some  in  flexion 
and  others,  simultaneously  in  flexion  and  extension.  Whence 
the  idea  of  using  the  following  explicit  approximations,  which 
should  be  valid  in  any  situation  : 
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where  the  coefficients  cf!^  and  c.  are  considered  as  constant 

i  J  t  .1 

throughout  the  redesign  phase. 


Because  it  is  no  longer  possible  to  select  appropriate  in¬ 
termediate  design  variables,  the  explicit  approximations  (33) 
cannot  be  obtained  by  merely  using  first  order  lav  I  or  series 
as  in  the  case  of  pure  bending  elements.  However  an  essential 
requirement  is  that  these  approximations  remain  correct  up  to 


J  . 


Che  first  order,  despite  the  fact  that  they  do  not  result  t row 
a  strict  linearization  process.  In  other  words  the  following 
equality  must  hold  : 
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This  condition  insures  that,  at  the  optimum,  the  solution  to 
the  explicit  approximate  problem  satisfies  the  (first  order) 
optimality  conditions  of  the  real  problem,  that  is,  the  appro¬ 
ximate  and  real  restraint  surfaces  have  the  same  tangent  plane 
(see  Fig.  1).  As  a  result,  the  GOC  approach  should  converge 
to  a  true  (at  least  local)  minimum  weight  design.  It  will  be 
shown  in  this  section  how  such  first  order  explicit  approxima¬ 
tions  can  be  obtained  for  various  types  of  behavior  constraints 
and  structural  models. 


4.1.  Displacement  Constraints 


The  key  idea  to  obtain  explicit  approximations  of  the  dis¬ 
placement  constraints  is  to  come  back  to  the  virtual  load  pro¬ 
cedure,  which  permits  decomposing  any  static  response  quantity 
into  the  contributions  of  each  element.  The  expression  (b)  can 
be  rewritten  in  the  more  general  form 
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Now,  for  a  rather  general  class  of  structural  models,  each  ele¬ 
ment  stiffness  matrix  can  be  assumed  to  have  the  following  ex¬ 
plicit  form  in  terms  of  its  design  variable  [  14  | 
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where  the  matrices  K  ^  P^  are  independent  of  the 
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are  assumed  to  be  constant  during 
dient  of  this  explicit  approximate 


P  c 


<P> 


aP  + 

l 


the  current 
constrain  t 


stage, 
i  s 


The  gra- 


(40) 


On  the  other  hand,  differentiating  (37)  and  inserting  the  result 
into  (36)  shows  that 
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Therefore  it  can  be  concluded  that  the  expressions  (38)  represent 
first  order  explicit  approximations,  in  that  they  restitute  the 
exact  values  of  the  constraints  and  their  first  partial  deriva¬ 
tives  at  the  design  point  a°  where  the  structural  analysis  is 
made  : 
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An  alternative  approach  is  to  employ  simple  linear  approxi¬ 
mations  in  the  reciprocal  design  variables  (9),  i.e.  the  first 
order  Taylor  series  expansions  (10).  Then  the  explicit  sub  pro¬ 
blem  exhibits  exactly  the  same  form  (14-16)  as  in  the  case  of 
thin  walled  structures  and  dual  methods  need  not  be  modified  as 
they  will  be  in  section  5.  Of  course,  as  previously  mentioned, 
the  convergence  of  the  overall  optimization  process  might  be 
lowered,  or  even  become  unstable,  in  the  case  where  the  bending 
behavior  is  dominant  in  most  of  the  elements.  The  reader  is  re¬ 
ferred  to  sections  6.1  and  6.5  for  numerical  examples  comparing 
the  explicit  approximations  (8),  (25)  and  (38). 

4.2.  Stress  Constraints 


The  situation  is  much  more  delicate  for  stress  constraints 
than  for  disp1 acement  constraints.  In  contrast  to  the  case  of 
thin  walled  structures  modeled  by  bar  and  membrane  elements,  the 
stress  matrices  are  no  longer  constant  in  the  case  of  flexion- 
extension  elements.  For  illustration  consider  again  a  flat  shell 
element.  The  stresses  in  the  upper  sheet  can  be  computed  in 
terms  of  the  generalized  displacements  by 
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where  o,  ,  are  matrix  representations  ol  the  total 

stress,  membrane  stress  and  bending  stress  tensors  in  member  k, 
and  T^,  are  the  corresponding  stress  matrices,  which 

are  independent  of  the  design  variables.  Each  displacement  com¬ 
ponent  admits  a  first  order  explicit  approximation  of  the  form 
(see  Eq.  33) 
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Therefore,  a  natural  choice  for  the  approximation  of  a  stress 
component  could  be  as  follows  : 
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The  expression  (43)  is  still  a  first  order  approximation. 

It  is  no  longer  separable,  but  keeps  a  simple  algebraic  form 
suitable  to  specialized  algorithms  such  as  those  based  on  dual 
methods.  However,  when  considering  the  true  stress  limitation 
(e.g.  upper  limit  on  the  Von  Mises  stress),  the  explicit  cons¬ 
traint  becomes  much  more  complicate.  For  these  reasons,  and 
before  finding  something  better,  it  has  been  decided  to  follow 
the  alternative  strategy  proposed  at  the  end  of  section  4.1, 
that  is,  to  employ  simple  linear  approximations  of  the  form  (10). 
To  compute  the  gradient  of  the  stress  constraint  with  respect  to 
the  reciprocal  variables  (9),  we  note  that,  from  the  equivalent 
Von  Mises  stress 
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2  x .  x  y  h  x .  y  x  i  x . 

l  i  J  j 
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I  12 

k 

+  6  T 


(46) 
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xy 


xy  3xi  ‘  k 


1  ,  (4 7  ) 


On  the  other  hand,  from  (44)  it  follows  that  tin-  derivative  <> 

as 

(f  ) 


any  stress  component  o  or  i  lias  the  f  o  t  in 


3  o ,  u 

_ x  =  T  ^9  _  s 

9  x .  k  3  x .  Lik  2 
t  1  x . 

l 

where  6^  is  the  Kronecker  symbol. 
4.3.  Frequency  Constraints 


( 4  H  ) 


Constraints  on  natural  frequencies  usually  consist  in  imposing 


1.  .  (  a  ) 


w >  0 

-J 


m 


1  owe  r  limits 


'ii?  (a) 


J 


(4'0 


f 


1 

I  9  . 


They  are  directly  written  in  terms  of  the  squares  of  the  fre¬ 
quencies,  because  these  quantities  naturally  appear  in  the 
eigenproblem  characterizing  the  structural  modal  analysis  : 


Kq.  -  wt  M  q  .  =0 
J  J  J 


(  ><>) 


in  this  equation  K  and  M  represent  the  stiffness  and  mass  matri¬ 
ces,  and  (q.,  j  =  1  , m)  are  the  modal  displace  me nts,  i . e .  ,  the 

J  .  2 

eigenvectors  solution  of  (50),  associated  with  eigenvalues  ,5. 
The  structural  mass  matrix  has  a  linear  form  in  terms  of  the 
design  variables  : 


n  n 

M  =  .  E  ,  M .  +  M  =  .  E  ,  a .  M .  +  M 
1=1  1  c  1=1  1  1  c 


(5  1) 


where  and  are  independent  of  the  design  variables.  M ^ 

denotes  the  mass  matrix  of  the  i th  element  when  a.  =  1.  M  re- 

i  c 

presents  the  contribution  of  the  non-s t ructu ral  masses,  such  as 
equipments,  fuel,  etc...  It  is  well  known  that  the  first  deri¬ 
vatives  of  the  frequencies  with  respect  to  the  design  variables 
are  given  by  |  see  for  example  Ref.  [  1  ],  section  4.2  ]  : 


3  oj '.  .  ...  9  K  . 

_J.  =  L_  q^(_i 

3a.  m .  j  J  a . 

l  j  J  l 


2  3Mi 
111  j  3a. 


)qi 


(52) 


where  m . 

J 


is  the  generalized  mass  of  the  j th  mode  : 


T 

m  .  =  q  .  M  q  . 

J  J  J 


( 5  5 ) 


The  way  to  derive  first  order  explicit  approximations  ol  the 
frequency  constraints  is  less  apparent  than  for  displacement  cons¬ 
traints.  In  this  report,  guided  by  the  work  done  in  Ref  (  2  | 
(section  4.2),  the  following  decomposition  of  the  eigenvalues  in 
terms  of  the  stiffness  and  mass  contributions  of  each  element  will 
be  used  : 


2 
u  ■ 
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iii 


q!  (K.  -  m2  M .  ) q  . 
J  i  J  i  J 


(54) 


=  m 
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whe  re 


m  . 
J 


(55) 


represents  the  contribution  of  the  non  structural  masses  to  the 


generalized  mass  m  j  (see  Lq.  hi).  By  taking,  account  of  the  ex¬ 
plicit  definitions  (37)  o l  the  stiffness  matrices  K  ^  and  (hi) 
of  the  mass  matrices  M.  ,  the  i  i  i  ;;  h  quality  explicit  approximations 
of  the  frequency  constraints  take  the  form  (>!)),  with 


o  m  j 

u.  =  ^t(l  +  -J-) 
.1  J  raj 


c<p) 

i  J 


T  M . 

J  t 


lp)qj 


P  = 


where  =  1  only  if  p  =  1  olid  is  0  otherwise. 


The  coefficients  c 


(P) 


and  the  modified  limits  u.  are 


i  J  J 

their  values  at  the  current  design  point.  lust  as  fo 

placement  constraints,  it  is  easily  verified  that  the 

(  3  K  )  remain  first  order  explicit  approximations  satis 

43).  As  a  matter  of  fact,  they  can  be  interpreted  as 

Taylor  series  expansions  in  terms  of  )  /  x  .  ,  )  /  x  T  ,  1  /  x  .* 

as  independent  variables. 


( r>  f>  ) 


1,2,3  ( 5  1  ) 

frozen  to 
r  the  d  i  s  - 
express)  mis 
f  y  i  n  g  (42, 
first  order 
considered 


4.4.  Buckling  Constraints 


Just  as  the  natural  frequencies,  the  critical  load  factors 
Xj  are  defined  through  an  e i ge np rob  1  e m 


K  q  .  -  '  -  S  q  .  -  0 

J  J  1 


(  r>  8  ) 


where  S  represents  the  geometric  stilfness  matrix  and  ( q - , 
j  =  I, m)  denote  tiie  eigenvector.-,  solution  ot  problem  (58),  asso¬ 
ciated  with  eigenvalues  a  .  I  he  pli  v  s  t  e  a  1  meaning  of  the  q  -  1  s 
is  that  of  displacements  in  the  jth  buckling  node  ,  for  a  criti¬ 
cal  load  factor  file  buckling  constraints  consist  in  imposing 

lower  limits  on  the  buckling  loads  : 


\  .  >  >  . 

.1  I 


j  =  I  ,  m 


<  )  4  ) 


In  this  work,  the  following  form  ol  th 
adopted  : 


taints  will  b i 


hL 


j 


2  I  . 

>  /  N  1  1  . 

hj(a)  —  -  ---  ( a  y  '  U  j  =  I,  m  (f.U ) 

because  it.  has  been  found  that  better  explicit  approximations 

are  generated  when  expanding  the  reciprocal  of  the  buckling 

loads  rather  then  the  A.'s  themselves. 

J 

The  stiffness  matrix  K  has  the  form  (23).  The  geometric 
stiffness  matrix  is  related  to  the  initial  stress  state  in  the 
elements  and  therefore  it  depends  implicitly  on  all  the  design 
variables  : 

n 

S  =  L I ,  S £  ( a )  (61) 

It  is  worth  recalling  that  the  matrices  S  ^  are  independent  of 
the  design  variables  for  a  statically  determinate  structure. 


As  explained  in  section  5.2  of  Ref  [ 
tives  of  the  buckling  loads  are  given  by 


the  tirst  d  e  r  i  v  a  - 


J  \  . 

_ l  = 

3  a  . 


1 


q  .  S  q  . 
J  J 


T  /Ki  _  ,  3S 
q  j  3a.  j  3a. 


)q  ; 


(62) 


In  opposition  with  the  static  and  dynamic  cases  previously  dis¬ 
cussed,  the  derivatives  appearing  in  (62)  are  not  directly  a  - 
vai table,  because  the  elements  of  the  geometrical  stiffness  ma¬ 
trix  are  function  of  the  stresses  acting  in  the  prebukling  state 


However,  by  assuming  that  the  terms  - 


3  S 
3a; 


are  negligible,  the  gra¬ 


dients  (02)  become  easily  computable.  This  assumption,  which  is 
typical  of  optimality  criteria  approaches  for  static  constraints 
amounts  to  not  taking  into  account  the  effects  of  structural 
redundancy  : 


3  S 

3  a . 

i 


=  0 


i  =  1 


(63) 


In  this  re port, the  following  decomposition  of  the  reciprocal 
buckling  loads  in  terms  of  the  contribution  of  each  element  will 
be  used  : 


-  q  .  K  .  q 

s  .  1  =  1  l  1  1 

1  1 


whu  re 


s  •  -  q  .  S  q 

I  .1 


Subs  ti Luting  the  explicit  del  ini  t  i  en  i  )  ;  tin-  s  r  i  !'  t  ness  n.a  - 

trices  into  (64  )  yields  tlu-  h  i  **  ti  quality  p  I  i  .  it  a  pp  ro  x  i  uni  - 

tions  of  the  buckling  load  constraints  in;.)  under  the  loin  (,3k) 

with 


,  ,  q*  Kf[,)  q  .  , 

cfP)=  - 1  a-l> 

j.  J  2  t 

s  . 

J  J 

Again,  the  coefficients  and  the  modi  tied  limits  u  .  are 

^  J  J 

frozen  to  their  values  at  the  current  design  point.  lust  as 
for  the  displacement  and  frequency  eons  t  r  a  i  n  l  s  ,  it  is  easily  ve¬ 
rified  that  the  expressions  (  18  )  am  i  i  r  s  t  .  r  .1  e  :  exp  I  i  e  i  t  appro- 
xi  mat  ions  of  the  constraints  (  nil )  ,  s  a  t  i  .  t  ■  i  n  t  • .  ,-•*),  provided 
that  the  assumption  (hi)  can  be  i  n  t  r  o  u  n  e  d  into  the  gradient 
(62)  (statical  de  L  e  rin  i  n  a  cy  )  . 


The  reason  why  the  buckling  e  .. :  t.  t  a  i  u  t  s  i  re  treated  by  ex¬ 


pand  i  ng  tlie  reciprocal  critical  loads 


1  q  .  6  4  )  is  that,  to r 


a  thin  walled  structure,  the  a  p  p  r  o  x  i  r.i  a  t  i  i  r  v  16  I  takes  the  I  or  in 

(  8  )  ,  where  the  coefficients  c  .  .  r  e  n  a  i  n  n  nsu.it  along  the  s  e  a  1  i  n  , 

1  ! 

line.  This  property,  which  is  well  detailed  in  section  3  o  t 
Ret  I  2  |  ,  has  been  found  very  important  for  the  case  <>  f  stress 

and  displacement  constraints.  C.'eume  t  r  t  o  a  1  1  v  it  means  that  t  lu¬ 
re  a  1  restraint  surtace  it  .  ( a )  =  ()  is  replaced  bv  a  tangent  sur- 
face  h  j  ( a )  =  0  at  its  point  of  intersection  with  tlu-  scaling, 

line  (see  F i g .  3  of  Kef.  |  2  |  )  . 


5.  SOLUTION  OF  T1IF.  EXPLICIT  PROBLEM 


From  the  foregoing  developments,  it  appears  t  li  a  t  ,  at  eaeli 
stage  o  t  Lite  optimization  process,  the  following  s  ubp  roll  I  cm  mu:.i 
be  solved 


minimize 
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l  l 

(68) 

subject  to 
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(69) 
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1  ph 
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^  a  . 
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(70) 

Unless  the  summation  on  p  in  (69)  is  restricted  to  a  single 
value  of  p  (i.e.  approximation  of  the  form  (25)  with  p  =  1,  2 
or  3),  it  is  no  longer  possible  to  find  intermediate  variables 
| i.e.  x^  given  in  (24)  1  in  terms  of  which  the  explicit  cons¬ 
traints  (69)  would  be  linear.  Therefore,  the  primal  solution 
of  problem  (68-70)  is  more  difficult  to  achieve  if  a  gradient 
projection  type  of  algorithm  is  employed  as  in  the  mixed  method 
developed  in  a  previous  work  for  thin  walled  structures  l  l  ,  1  !>  |. 

The  expressions  (69)  are  still  explicit  and  they  continue  to  ex¬ 
hibit  a  simple  algebraic  form.  Consequently  a  general  purpose  op- 
t.  i  m  i  z  a  t  i  o  n  a  l  go  r  i  Lhm  s  uch  as  N  KWS II MT  |  16  |  could  e.isi  lv  In'  ad.iptnl  to 
take  the  constraints  (69)  into  account.  However,  because  they 
are.  still  additively  separable,  resorting  to  dual  methods  remains 
probably  the  best  strategy,  just  as  in  the  case  ol  thin  walled 
structures  [1,  2,  6,  8  ]. 

The  minimization  problem  (68-70)  can  be  solved  efficiently 
as  an  auxiliary  maximization  problem  in  the  m  1  a  g  r  a n  g i a n  multi¬ 
pliers  r j  associated  with  the  explicit  behavior  constraints  (09). 
This  dual  problem  reads  as  follow  |  8  ]  : 

n  m 


maximize  H  (  r )  =  .  L  .  C. .  a.(r)  +  .5'  r.  g  .  ( r )  (71) 
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(7  1) 


The  primal  variables  .  t  .  (  r  )  are  related  t  a  t  1 ;  <  •  dual  variables  r  . 

L  I 

through  the  loll  viw  i  a-,  ■  j  u  a  s  i  -  u  u  e  s  l  r  a  i  in  d  mi  nimi  /.atimi  p  r  <>  l>  I  i  ■  m 

(see  foi  example  Ke  I  jo  i  )  : 
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u j  )  (  74  ) 


subject  to  a. 


a  .  a  . 

l  i 


{  7  1)  ) 


Because  of  the  separability  ol  this  problem,  it  can  be  decompose 
into  u  one-dime  n  siouai  ra  i  n  i  mi  ?.  a  t  1  yn  prop  le  ms  o  j  the  form 
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~i 


a  .  a  . 
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(7b) 


where 


(  p  )  al 

C  : 1  =  p  ■  -  ,  r  .  c  : 

i  J  =  1  J  i.l 


(P) 


(77) 


Setting  to  zero  the  first  derivative  o  1  t  h  ••  single  variable 
function  appearing  in  |:7b),  it  is  .ecu  that  the  n  "active"  de¬ 
sign  variables  can  be.  obtained  b  solving  the  nonlinear  alge¬ 
braic  equations 


:>  ‘i 


no 


p  = i  a?’  i 
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I  ,  n 


(78) 


Note  however  that  when  solving  ,  be  .-ide  constraint  (7b) 

must  be  taken  into  account  and  'bat,  if  several  values  satisfy 
equation  (78),  the  one  that  mi  n  i  mi  /.i.  e  ti.o  function  (7b)  must  be 
retained.  Standard  techniques  can  bo  employed  to  deal  with  such 
a  simple  one  - d i me  ns i ona i  problem  ( s c e  t  . g .  |  1  J  |  )  . 


In  many  cases,  as  previously  men  l  i  ri’.t  d  ,  at  least  one  ol  the 
terms  in  the  summation  ■•!!  p  in  (  •'  i>  >  .i  n.i  (78)  disappears,  and 
it  is  possible  to  treat  the  problem  analytically.  As  an  illus¬ 
tration,  consider  again  the  case  o I  a  I  1  a  t  shell  element,  where 
the  term  in  p  -  /  is  missing, 
form  (omitting  the  index  i)  : 


i‘  rub  1  e  m 


(7b)  exhibits  then  the 
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whi  le  equation  (it.)  is  r  »  [  low 

,(;i  .  t  . ) 
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or,  setting 
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(82) 


The  constant  K  is  known  to  be  t  i  v  e  ,  and  e  v*  rv  thing  depends 

thus  on  tlie  sign  of  tin-  con, cants  ^  ^  and  0  ^  \  which  repre¬ 
sent  the  contribution  ot  the  membrane  and  the  bending,  respec¬ 
tively.  In  general,  equation  (82)  admits  two  solutions,  and  the 
one  that  renders  mini  m  a  i  the  (unction  (,/“)  must  he  adopted. 


From  the  foregoing,  developments,  it  appears  that  the  dual 
function  (71)  can  be  c  ■  ■  tin  i  d  e  l  ••  d  as  i  function  of  the  dual  vari¬ 
ables  only.  The.  dual  problem  stated  in  (71,  72)  exhibits  an 
attractive  feature,  namely  ,  'll  is  ,i  quns  i  -  uncons  t  rained  problem, 
because  taking  care  of  fin  no,,  n  e  g  a  t  i  v  i  f  v  constraints  (72)  on  the 
dual  variables  is  *  t  to.  i  -tit  i  \  •-  •  .!  .  i :  i  .  ■  u  ..••me  non nog a  tivc  dual 

variables,  the  c  o  r  i  e  s  p  o  mi  .  n  c  pi  i  ,i-i  i  v .  1 1  i  i  it  1  e  are  compu  ted  by  sol¬ 
ving  (7b)  or  (78>  and  the  :» r  i  m..  j  is.  I  Mint:,  are  evaluated  by  u- 
sing  ( 7  J  )  .  The  dual  Ciinc  t.  inn  i  -  ■  i  .ci.l  its  gradient  (78)  are  then 
directly  known  and  a  I  c  a  :•  i  h  i  e  a  ..  .•  o  n  t  d  i  recti  on  can  therefore  be 
determined.  In  a  seccinl  ordci  algorithm,  the  hessian  matrix  of  the 


dual  function  has  to  hi-  t,  oni  puti-d 

1 8 


1  i  ' 1  k 


From  (71),  it  c  o  me 


n  ') 
i  i  p 


i 


i  a  . 

t 


(89  ) 


and  by  differentiating  the  definition  (78)  of  a.(r)  with  respect 

to  r .  ,  it  follows  that  . 
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Oil  tlie  other  hand,  i  l  is  oas  1  1  v  sac:,  f  ruin  (77)  that 


Regrouping  terms,  it 
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where  the  summation  or:  i  is  restricted  to  the  n  active  primal 

variables,  that  is  ,  the  pi.  that  arc  not  lixed  to  a  lower  or  an 

j  a  • 

upper  limit  (of  course  - - =  0  lor  a  passive  variable).  Knowing 

‘ 1  k 

the  gradient  (73)  and  tiie  hessian  matrix  (87)  furnishes  the 
Newton  searcli  direction 


=  -  H 


The  next  dual  point  is  then 


given  hv 


where  a  is  the  step  1  .•  ng  -J: 
often  a  regular  Newton  t 
lected,  however,  the  value 


;  n  e  jirm  tion  /  .  Most 
i  i  in  ( 8  9  )  |  is  sc 

■•••<  I  i  e  s  be  lowered  to 


prevent  one  of  the  dual  variable  ..  t  rom  he  coming  negative 


The  second  order  dual  optimist  t  i  tap  i  c  me  n  t  e  d  in  SAMO.KF  |  IK  | 
has  been  esiecial ly  devised  so  that  it  socks  the  maximum  ot  the 
dual  function  bv  operating  in  a  s  e  i|  ue  n  i  <■  of  dual  s  u  b  s  p  a  c  e  s  w  i  t  h 
gradually  increasing  d  i  mo  psii.-n.  In  this  way,  the  effective  di¬ 
mensionality  o  1  the  in ,'i  \  i  n,  i  /it  inn  problem,  never  exceeds  the  num¬ 
ber  of  active  behavior  on?  I.  i  ,i  i  n  ts  ,  which  correspond  to  non¬ 
zero  dual  variables.  Past  experience  with  thin  walled  structures 
indicate  that  this  number  is  relatively  small  in  practice,  w  h  i  c  ti 
explains  the  remarkable  el  I  i  >  i  emv  of  the  dual  method  |  6  j  or 
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variables  in  terms  o  1  t  in  d  u  a  I  m,-  i  jHis  |  sit  K  <| ..  (  7  6  ,  7 8 )  ]  and 

the  formula  employed  I  ...  r  t  in  ii  .•  ••  .  a  i.  ni-iiiix  |  flee  Lq  (  H  7  )  1  .  It 
should  however  be  recognized  t  ha:  ..  u e  i  i  e  a  1  difficulties  might 
occur  due  to  the  lack  o  t  co ve  m  t  y  t  ti:  e  explicit  constraints 
(69).  Indeed  it  is  possibl.  that  :  Ik-  solution  u  I  l.q  (76)  lor  a 
primal  variable  suddenly  ju.ips  !  to,.:  a  l  ixed  lower  or  upper  bound 
to  a  tree  value  when  the  dual  variables  are  slightly  modified. 
This  phenomenon,  which  means  that  the  dual  function  first  deri¬ 
vatives  are  discontinuous  alone,  sc. me  surfaces,  needs  more  inves¬ 
tigation. 


6.  NUMERICAL  APPLICATIONS 


All  the  examples  presented  involve  rather  sophisticated 
flat  shell  elements  that  are  characterized  by  a  displacement 
field  cubic  in  extension  and  quintic  in  flexion  (hybrid  qua¬ 
drangular  flat  shell)  [  18  ]  . 

6.1.  Cantilever  Beam  with  End  Moment 


The  first  example  is  concerned  with  a  cantilever  beam  loaded 
with  a  concentrated  moment  at  its  free  end  1  see  Fig.  2  ] .  This 
problem  has  been  previously  solved  numerically  by  PRASAD  and 
HAFTKA  |  19  ]  using  an  extended  interior  penalty  function  for¬ 
mulation  for  a  beam  having  the  following  properties  :  length  = 

10  in,  width  =  1  in,  applied  moment  =  540  in. lb  ,  Young's  modu¬ 
lus  =  107  psi,  Poisson's  ratio  =  0.3  and  mass  density  =  0.3  lb/in 
There  exists  also  an  analytical  solution,  which  was  obtained  by 
HAUG  |  20  ]  .  The  analytical  optimum  design  for  a  displacement 
limit  of  0.5  in  at  the  free  end  and  specified  allowable  stresses 
of  30,000  psi  is  as  follows  (see  Fig.  2)  : 


a ( x )  =  0.30 

for 

0  <  x  <  2 .  3 

(90) 

a(x)  =  0.244  ^ 

for 

2. 3  <  x  <  10 

(9  1  ) 

It  is  worth  noticing  that  the  beam  being  statically  determinate, 
the  redesign  relation  (31)  is  exact  (with  p  =  3)  and  reduces  to 
(90,  91).  Therefore  the  optimal  solution  should  be  generated 

after  one  structural  analysis  only. 

The  structure  is  discretized  using  20  quadrilateral  plate 
bending  finite  elements  as  indicated  in  Fig.  2.  There  is  no 
design  variable  linking  and  the  problem  involves  thus  20  inde¬ 
pendent  design  variables.  The  problem  was  first  solved  by  using 
the  cubic  expansions  (25)  with  p  =  3  and  the  dual  algorithm  des¬ 
cribed  in  section  5.  As  previously  stated,  because  the  structure 
is  statically  determinate,  the  explicit  problem  (68-70)  is  exact 
and  the  optimum  design  is  obtained  in  one  single  analysis  what¬ 
ever  may  be  the  initial  thickness  a°.  Then,  the  same  example 


i _ 


h  )  will'  a  ii  i  n  i  - 


was  again  solved  bv  using  the  linear  ex  pans i on 

.  .  <  I  .  .  .  . 

U  .1  l  thickness  a  0  .  >  in  ,  r  e  s  u  I  t  i  u  g  in  an  appro  x  j  m  ate  e  x  p  J  in  t 

problem  of  the  form  (14-lb).  The  iteration  history  is  given  in 

Table  1.  It  can  be  seen  that  ,  s  in  p  r  i  s  i  ng  L  y  enough,  the  use  of 
simple  linear  expansions  gives  i  i  .  e  to  satisfactory  results. 

The  two  final  designs  obtained  are  listed  in  Table  2,  toge¬ 
ther  with  the  analytical  solution  corresponding  to  Kqs.  (90,  91) 

and  the  numerical  results  ot  Re  I  |  19  j  .  A  close  agreement  bet¬ 
ween  all  the  final  opt  mu!  des  i  i  an  lie  observed. 

b  .  2  .  S imply  Supported _ S  q  uu  re  I'uii'  u  i  t  ii  Deflection  Constraint 

The  second  example  consist  2.  ii  it.  i  n  i  m  i  z  i  n  g  the  weight  of  the 
simply  supported  square  plate  shown  <>n  !■  i  g .  3.  It  is  subjected 

to  a  concentrated  load  of  1 000  ki  »  t  its  center,  where  the  de¬ 
flection  is  limited  to  0.0.'  m .  >\  :•  visi.ie  try  only  a  quarter  of 

the  plate  has  to  be  analysed.  i he  mesh  involves  29  plate  ele¬ 
ments  and  179  degrees  of  freedom.  ihe  material  properties  are 
as  follows  :  Young's  modulus  f  .  .Jo*1  N  in'  ,  Poisson's  ratio 
v  =  0.3  and  weight  density  V. .  :•  * .  Hie  dimensions  of 


the  plate  are  10  m  b;  0 

m . 

i  ii  <. 

.  .  o 

.*  tiiiU.il  tli  i  -  kiu  ss  is  a 

=  0.0  9  m 

and  the  minimum  thickness 
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roll  Jem  in¬ 

volves  only  one  behavior 
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jurious  (  f  1  ) 
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is  illustrated  in  Fig.  5. 

I 

t  i  - 

v-  tfh  ;■  ;  i 

n  t.  i  n  g  o  u  t 

t  h  a  t 

only 

)0  reanalyzes  are  s u f  f i c i 

f  Ii  t 

f  , 

>,  •  ii  e  r  a  t  e  a 

nearly  converged  so- 

lution,  while  in  a  si  mi  la 

r  P 

r  1 1  b  1 

em  solved 

in  Ref  |  12 

1  by 

usi  n  y. 

p  =  2  in  the  redesign  relations  (31),  almost  90  iterations  were 
r  eq  u  i  re  d  . 


b.J.  Simply  Supported  Square  I’l.itc  with  Frequency  Constraint 


Attention  is 
shown  in  Fig.  4  . 


now  focused  on  the  simply  supported  square  plate 
This  example  is  taken  from  Kef  |  2  I  |  and  it.  is 


I 

■  4 


30  . 

concerned  with  the  minimum  weight  design  of  the  plate  subject 
to  a  natural  frequency  constraint.  The  dimensions  of  the  plate 
are  10  in  by  10  in  and  its  material  properties  are  Young's  mo¬ 
dulus  =  30  x  1  0  ^  lb/in^  and  Poisson's  ratio  =  0.3.  In  Ref  I  31  | 
i't  is  stated  that  the  uniform  plate  with  pco  =  1400  is  taken  as 
the  initial  estimate  to  the  optimization  problem,  where  p  is  the 
mass  density  and  w,the  minimum  fundamental  frequency.  Therefore, 

in  the  present  study,  assuming  steel  material  with  miss  density 
3 

p  =  0.283  lb/in  ,  the  minimum  frequency  was  chosen  as  =  id/  2  ;i  - 
11.2  Hz.  Iteration  history  data  are  presented  in  Table  4  for 
three  cases  differing  by  the  minimum  thickness  constraint  (a  =  0. 1 
0.05  and  0.001  in).  The  initial  design  in  cases  1  and  3  corres¬ 
ponds  to  a  uniform  thickness  a°  =  0.2  in,  while  in  case  2,  n° 
equals  0.12  in.  Final  uesignsare  illustrated  in  Fig.  4.  By 
symmetry  only  a  quadrant  of  the  plate  has  to  be  analyzed  and 
designed.  It  can  be  observed  that  the  design  obtained  in  Kef 
121  ]  is  different  from  the  design  achieved  in  this  study  (case  I) 
However  both  designs  have  about  the  same  weight  (0.765  lb  and 
0.752  lb).  It  is  also  worth  pointing  out  that  the  design  of  pla¬ 
tes  in  bending  is  an  analytically  difficult  problem,  leading  to 
multiple  local  optima  and  suggesting  that  the  optimal  plate  should 
be  made  up  of  an  infinite  number  of  stiffeners  |  22  )  . 

6.4.  I-beam  Structure 

The  third  example  involves  the  I-beam  structure  schematized 
in  Fig.  5.  The  problem  consists  in  minimizing  the  weight  ol  the 
beam  while  imposing  lower  bounds  on  the  frequencies  of  the  three 
first  eigenmodes  :  flange  flexion,  torsion  and  web  flexion.  De¬ 
tailed  data  can  be  found  in  Ref.  (  2  ] .  In  a  first  optimization 
exercise,  a  pure  membrane  model  was  employed.  It  involves  35 
second  degree  displacement  elements,  including  10  fictitious 
diaphragms  (without  masses).  These  dummy  members  are  introduced 
to  obtain  a  satisfactory  representation  of  the  torsional  mode. 

Only  5  analyses  are  sufficient  to  generate  an  optimum  design  for 
this  membrane  model. 


3  I  . 

However,  when  this  final  design  was  analyzed  by  using  a  more 
accurate  model  made  up  of  flat  shell  elements,  the  torsional  fre¬ 
quency  (inode  2)  was  seen  to  be  violated  by  10  7.  Therefore  tin 
problem  was  again  solved  with  this  n  e  w  ru  <  •  d  <  I  ,  by  • <  s  o  r  t  i  m*  ' 1 '  the 
theory  proposed  in  section  4.3  and  to  the  dual  optimizer  descri¬ 
bed  in  section  5.  Iteration  history  data  are  illustrated  in 
Fig.  5  and  the  final  designs  are  given  in  Table  3  for  both  finite 
element  models  of  the  I-beam.  It  can  be  seen  that  the  use  of 
flat  shell  elements,  although  yielding  slower  convergence,  gives 
satisfactory  results. 

6.5.  U-beam  Structure 

In  an  attempt  to  consider  a  case  where  both  flexion  and  ex¬ 
tension  loadings  play  an  equally  important  role,  the  U-beam 
structure  depicted  in  Fig.  6  was  optimized.  The  2000  kg  toad 
acting  at  the  tipproduces  torsion  of  the  beam.  As  a  result  the 
upper  flange  behaves  mainly  in  flexion,  the  web,  in  extension, 
and  the  lower  flange,  both  in  flexion  and  extension.  For  sim¬ 
plicity,  only  displacement  constraints  are  considered  :  the  tip 
deflection  is  limited  to  ().!  m,  while  the  relative  lateral  dis¬ 
placement  is  limited  to  0.005  m  (see  Fig.  6).  Three  different 
methods  are  employed,  which  ditfei  by  the  explicit  approximations 
used  for  tlie  displacement  constraints  : 

use  I  :  full  expansion  (  3  ,  yielding  problem  (68-70  )  with 

p  =  1  and  3  only  i,  tin  term  in  p  =  2  being  zero) 

case  2  :  linear  expansion  2  >  witii  p  -  !),  yielding  problem 

(14-16) 

case  3  :  cubic  expansion  (25  with  p  =  3 )  . 

1  he  iteration  history  data  are  given  in  Table  6  and  the  final 
designs  in  Table  7. 

It  is  interesting  to  notice  that  the  three  methods  give 
similarly  good  results,  wh Leh  suggests  that  the  use  of  simple 


linear  Taylor  series  expansion  in  terms  of  the  reciprocal  di-si  j;n 
variables  might  be  a  good  strategy  for  most  of  the  structural 
optimization  problems.  This  idea  was  previously  stated  by  AUSTIN 
[  23  ]  ,  but  not  believed  valid  by  many  people  when  he  published 
his  paper  ! 


!  i  . 
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it  is  now  widely  recognized  that  a  powerful  and  rather  ge¬ 
neral  approach  to  structural  o ;  >  t  i  in  i  z  a  i  i  •  >  n  is  achieved  bv  repla¬ 
cing  the  original  problem  with  a  sequence  of  explicit  approxi¬ 
mate  problems.  This  approach  was  initially  conceived  for  thin 
walled  structures  modelled  bv  bar  and  membrane  finite  elements, 
as  well  in  the  context  of  optimality  criteria  techniques  as  in 
the  framework  of  mathematical  programming  methods  using  approxi¬ 
mation  concepts.  It  has  been  extended  in  this  work  to  deal  with 
structural  systems  made  up  of  beam,  plate  and  flat  shell  element!., 
with  behavior  constraints  placed  on  displacements,  stresses,  na¬ 
tural  frequencies  and  critical  buckling  loads.  The  method  p re¬ 
sented  uses  a  second  order  dual  algorithm  to  solve  each  explicit 
subproblem.  The  convergence  properties  are  independent  of  the 
number  of  design  variables,  which  is  typical  of  optimality  cri¬ 
teria  types  of  approach  as  well  as  linearization  techniques  in 
mathematical  programming.  As  a  result  large  structural  systems 
can  be  treated  at  the  expense  of  a  few  finite  element  analyses. 

Several  difficult  points  remain  to  be  clarified.  First  it 
is  not  well  known  whether  the  explicit  problem  can  be  solved 
conveniently  in  any  case,  because  its  lack  of  convexity  might 
lead  to  discontinuity  in  the  dual  function  gradient.  Secondly, 
the  use  of  simpLe  linear  expansions  with  respect  to  the  recipro¬ 
cals  of  the  element  transverse  sizes,  gives  rise  to  satisfactory 
results  in  many  cases.  This  suggests  that  more  complicate  expli¬ 
cit  approximations  such  as  those  proposed  in  this  work  mi  g  lit  not 
be  necessary.  Finally,  there  remains  the  question  of  the  stress 
constraints,  for  which  it  is  difficult  to  generate  adequate  first 
order  explicit  approximations.  Also  s t r e s s - r a t i o i ng  algorithms 
should  be  considered  lor  extension- flexion  elements,  because  they 
correspond  to  much  more  inexpensive  zero  order  approximations. 
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Table  2  Final  Designs  for  Cantilever  Beam 


E  lenient 
No. 


Analytical 
Ref  [  20] 


Thick  ness  (in) 


Ref  [  19  ] 


Linear 
Expans  ion 


ub  i  c 
Expansion 


1 

2 

3 

4 

5 

6 

7 

8 
9 

1  0 
1  1 
1  2 
1  3 
14 
1  3 
1  6 
1  7 
1  8 
1  9 
20 


0 . 3 
0.  3 
0. 3 
0 . 3 
0 . 3 

0.3142 
0.3276 
0. 3396 
0.3503 
0.3602 
0.3693 
0.3778 
0 .3858 
0.3933 
0 . 4004 
0.407  1 
0.4135 
0.4197 
0.4255 
0.43  12 


0 . 3004 
0 . 3006 
0.3018 
0 . 3056 
0.3109 
0.3167 
0.3231 
0 . 3298 
0 .3367 
0.3437 
0. 3508 
0.35  79 
0 .3652 
0 .3726 
0 .3803 
0.3881 
0 . 3960 
0.4040 
0.4119 
0.4197 


0  .  3 
0.  3 
0.  3 
0 . 3 
0.  3 

0.3167 
0.3287 
0.3327 
0 . 3350 
0.3379 
0. 3450 
0.3574 
0.3711 
0.3851 
0. 3989 
0.4121 
0.4246 
0.4358 
0.4446 
0.4533 


0 . 3 
0  .  3 
0 . 3 
0  .  3 
0  .  3 

0.3136 
0.3270 
0. 3389 
0.3497 
0.3596 
0.3687 
0.3771 
0.3851 
0.3926 
0.3996 
0.4063 
0.4126 
0.4182 
0.4225 
0 .4268 


Weight 

(lb) 


1  .0750 


1.0453 


1.0768 


1.0797 


No  •  of 
Ana  lyses 


/ 


8 


2 


Table  4  : 


Iteration  History  Data  foi  Square  Plate 
with  Frequency  Constraint 


Ana  lysis 

No. 

Case  1 

• 

a  =  0  .  1  in 

Case  2  : 

a  =  0.05  in 

Case  3  : 

a  =  0.001 

Weight 

(lb) 

Frequency 

(Hz) 

Weight 

(lb) 

Frequency 

(Hz) 

Weight 

(lb) 

Frequeni 

(Hz) 

1 

1.415 

19.6 

0.849 

11.8 

1.415 

i  y .  5 

2 

1.097 

16.0 

0 . 786 

11.5 

1  .097 

16.0 

3 

0.913 

13.6 

0. 750 

11.5 

0 . 902 

13.7 

4 

0.812 

12.2 

0. 708 

11.7 

0.778 

12.4 

5 

0.771 

11.5 

0.658 

11.6 

0.696 

11.9 

b 

0.757 

11.3 

0.63  1 

11.4 

0.647 

11.5 

7 

0 . 753 

11.2 

0.617 

11.3 

0.622 

11.4 

8 

0.752 

1  1  2 

0.605 

11.3 

0.604 

11.4 

9 

0.752 

11.2 

0.596 

11.3 

0.588 

11.4 

1  0 

0.589 

11.3 

0.5  74 

11.3 

1  1 

0.583 

11.3 

0.559 

11.4 

1  2 

0.577 

11.3 

0.542 

11.4 

1  3 

0.52  1 

11.3 

IBM  370-158 

CPU  time 

2  1  8 

276 

307 

(sec) 

_ 

nge  d 


Table  5  :  Final  Designs  for  I-beam  Structure 


thickness  (mm) 


membrane  model 
flat  shell  model 


15.80 

16.59 


Upper  Flange 

17.30  11.67  6.143 

15.09  10.08  5.090 


1.813 

1  .36  3 


5.101 

3 .602 

3.329 

3.294 

1  .  9  9  7 

2 .056 

4.173 

3.965 

3.936 

2.313 

S  upported 


15.40 

20.51 


17.01 

17.98 


Lower  Flange 

11.55 

12.21 


6.074 

6.480 


1.79  7 
1.91  0 


Table  6  :  Iteration  History  Data  for  U-Beam  Structure 


Table  7 


Final  Designs  tor  II -lie  am  Structure 


Upper 
F 1 ange 


Web 


Lower 
F  1  a  n  g  e 


Thickness  (mm) 


Case  1 

:  Full  Expansion 

Case  2 

:  Linear 

V  y. pans  i  n 

Case  i 

:  Cubic  E 

X  p  tins  i  on 

2  1.6  ! 

r  1 

16.6 

1 

10.9 

6.00 

1 

2.29 

2i.2  : 

i6.4  ;  • 

11.2  i 

6.67 

1 

2.27 

2  1.2  ' 

L  i 

_ f _ 

16.4 

1 

1 

10.8  1 

4.82 

1 

1 

_ l _ 

2  .  39 

t 

f  « 

r  ( 

r  i 

r 

1  4  .  a  , 

11.8  i 

8.  hi  , 

a  .  3  3 

t 

2  7.1 

14.0  !  . 

11.1  > 

• 

7  .  1  1  !  ' 

3.24  ' 

1 

2  7.0 

17.0  ' 

i3.8  : 

9.96’ 

6.00  1 

2  7.1 

_ _ J _ 

-  -- _ 1 _ 

J _ 

1 

** 

(24.6  ; 

’22.5  ! 

19.4 

T  — 

r23.6  ; 

17.9 

25.0  1 

23.4  :■ 

20 . 6 

1 

1  \ 

2  3.6  1 

17.9 

1  COMPARISON  OF  REAL  AND 
riL7.  1  APPROXIMATE  RESTRAINT  SURFACES 


finite  element  model 

;  i" 

it 

«* -  K)"-  -  * 

20  bending  plate  elements 
220  degrees  of  freedom 

FIG.  2  CANTILEVER  BEAM  WITH  END  MOMENT 


final  design(f/4)  thicknesses  (m) 


-  f/O- —  ' 

symmetrical 


SIMPLY  SUPPORTED  SQUARE  PLATE 
WITH  DEFLECTION  CONSTRAINT 


